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Abstract 


The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument on the Cloud- 
Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) spacecraft has provided 
global, high-resolution vertical profiles of aerosols and clouds since it became operational on 13 
June 2006. On 14 June 2006, the NASA Langley Research Center (LaRC) High Spectral 
Resolution Lidar (HSRL) was deployed aboard the NASA Langley B-200 aircraft for the first of 
a series of 86 underflights of the CALIPSO satellite to provide validation measurements for the 
CALIOP data products. To better assess the range of conditions under which CALIOP data 
products are produced, these validation flights were conducted under both daytime and nighttime 
lighting conditions, in multiple seasons, and over a large range of latitudes and aerosol and cloud 
conditions. This paper presents a quantitative assessment of the CALIOP 532 nm calibration 
(through the 532 nm total attenuated backscatter) using an internally calibrated airborne HSRL 
underflight data and is the most extensive study of CALIOP 532 nm calibration. Results show 
that average HSRL and CALIOP 532 nm total attenuated backscatter agree on average within 
2.7% ± 2.1% (CALIOP lower) at night and within 2.9 % ± 3.9% (CALIOP lower) during the 
day., demonstrating the accuracy of the CALIOP 532 nm calibration algorithms. Additionally, 
comparisons with HSRL show consistency of the CALIOP calibration before and after the laser 
switch in 2009 as well as improvements in the daytime version 3 calibration scheme compared 
with the version 2 calibration scheme. Potential systematic uncertainties in the methodology 
relevant to validating satellite lidar measurements with an airborne lidar system are discussed and 
found to be less than 3.7% for this validation effort with HSRL. Results from this study are also 
compared to those from prior assessments of CALIOP calibration and attenuated backscatter. 



1. Introduction 


The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite was 
launched in April 2006 (Winker et al., 2010). The main objective of the CALIPSO mission is to 
provide a global, multi-year data set of cloud and aerosol spatial and optical properties from 
which to assess uncertainties of aerosol direct and indirect effects on climate forcing and cloud- 
climate feedback (Winker et al., 2007, Winker et al., 2009). To address this objective, the 
primary payload on the CALIPSO satellite is a two-wavelength and polarization-sensitive elastic 
backscatter lidar, the Cloud- Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument. 
Since becoming operational in June 2006, CALIOP has provided unprecedented observations of 
aerosols and clouds in the atmosphere. This paper presents a quantitative assessment of the 
532 nm calibration by comparing CALIOP total attenuated backscatter profiles with coincident 
measurements acquired by the NASA Langley Research Center (LaRC) airborne High Spectral 
Resolution Lidar (HSRL). 

Since the launch of CALIPSO, several studies have been conducted validating the CALIOP data 
products. Some promising results were published by McGill et al. (2007) qualitatively 
comparing measurements from CALIOP to the Cloud Profiling Lidar (CPL) deployed on the 
ER-2 aircraft. McGill et al. examined the vertical distribution of clouds measured coincidently 
by the two instruments and assessed the minimum detectable backscatter from CALIOP. 
However, it concerned only a few profile comparisons showing quantitative agreement but does 
not provide a comprehensive quantitative assessment of CALIOP calibration accuracy. 


Furthermore, the CPL attenuated backscatter must be calibrated in the same manner as CALIOP 



while the internal calibration of the HSRL aerosol backscatter provides a truly independent 
calibration assessment. 

The airborne HSRL, developed by NASA LaRC (Hair et al., 2008), has been deployed on ten 
field experiments to-date, logging over 800 hours on 240 flights of the LaRC King Air B-200 
aircraft. Many of these ten field missions have included CALIOP validation flights. The HSRL 
technique (Piironen and Eloranta, 1994; Hair et al., 2001) allows the airborne HSRL used in this 
study to independently measure the aerosol backscatter and extinction profiles at 532 nm. The 
airborne HSRL also measures backscatter at 1064 nm using the standard backscatter technique 
and is polarization sensitive at both 532 and 1064 nm. The data products from the airborne 
HSRL are produced at high accuracy, with no assumptions in the calibration and few assumptions 
in the aerosol retrieval algorithms, and thereby offer an accurate and completely independent 
means by which to validate fundamental CALIOP products. Because of the importance of the 
532 nm total attenuated backscatter calibration, this paper focuses solely on this fundamental 
product. This paper also establishes a paradigm for validation of satellite lidars and other 
spaceborne instruments that provide nadir-only measurements: a systematic series of aircraft 
flights along the satellite track employing instruments with appropriate sampling geometries and 
measuring techniques to get an accurate and statistically significant database to validate the 
satellite products. HSRL validation of the CALIOP level 1 1064 nm attenuated backscatter and 
level 2 products are the subjects of future publications. 

2. HSRL Validation Flights 

The HSRL has completed 86 successful validation underflights of the CALIPSO satellite through 


2009, acquiring 116 flight hours of data along CALIPSO orbit tracks. The average time that 



HSRL spent along each CALIPSO track was 1.4 hrs + 0.8 hrs. In this time, the HSRL covers an 
average of 385 km along track per flight at an average ground speed of llOm/s. Figure 1 shows 
the locations of the CALIOP validation flight tracks flown with the airborne HSRL. These flights 
cover a wide seasonal and latitude range as well as a wide variety of aerosol and cloud 
conditions. 

As mentioned earlier, CALIOP validation flights were often flown by HSRL during field 
missions. Some of these missions were focused solely on CALIOP validation, however most of 
which had multiple objectives. Table 1 summarizes the specific mission information and further 
demonstrates the wide seasonal and spatial sample the HSRL validation flights encompass. The 
first such mission was the CALIPSO-CloudSat Validation Experiment (CC-VEX), where the 
HSRL was based out of NASA LaRC in Hampton, VA, and flew primarily along the US Eastern 
Seaboard. HSRL then participated in several other field studies that included CALIOP validation 
flights: the Texas Air Quality Study (TexAQS) - Gulf of Mexico Atmospheric Composition and 
Climate Study ( GoMACCS ) based out of Houston, TX (Parrish et al., 2009); the Cumulus 
Humilis Aerosol Processing Study (CHAPS) based out of Ponca City, OK (Berg et al., 2009); the 
CALIPSO and Twilight Zone (CATZ) field campaign based out of NASA LaRC; a mission 
dedicated to CALIOP validation flights based in the Caribbean islands; and the Arctic Research 
of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) spring and 
summer deployments based out of Barrow, AK and Yellowknife, NWT, Canada, respectively 
(Jacob et al., 2010). HSRL also conducted a special series of nighttime flights based out of 
Hampton, VA, to verify the long term stability of the CALIOP calibration. These flights covered 
the period when the CALIOP laser was intentionally transitioned from its primary laser to its 
backup laser due to a loss of pressure in the primarily laser (Hunt et al., 2009) and are discussed 



further in section 5.3. Most recently, HSRL participated in the Routine ARM Aerial Facility 
(AAF) Clouds with Low Liquid Water Depths (CLOWD) Optical Radiative Observations 
(RACORO) campaign based out of Ponca City, OK. In addition, flights not associated with a 
specific mission were occasionally conducted during transit flights to or from NASA LaRC and 
other destinations (denoted by “Other” in Table 1). 

3. Comparison Methodology 

3.1 CALIOP attenuated backscatter 

The CALIOP data products are divided into level 1 and level 2 categories as per the standard 
NASA Earth Observing System nomenclature.: level 1 and level 2 (King et al, 2004). The level 1 
products are geolocated and calibrated profiles of total attenuated backscatter coefficients at 532 
and 1064 nm and the perpendicular polarized component of the total attenuated backscatter 
coefficients at 532 nm. Total attenuated backscatter, denoted by J3'(r ) , is the sum of the parallel 
and perpendicular attenuated backscatter profiles as defined in equation 1. This paper deals 
solely with the 532 nm attenuated backscatter so no wavelength subscript is used on J3'(r ) . The 
parallel attenuated backscatter is defined as the product of the backscatter coefficient, /? N '(r) , and 

the two-way attenuation of the atmosphere, T 2 (r). This quantity is also written in terms of the 
measured parallel power, P\\(r), range, r, and non-range dependent parameters (the system 
calibration constant), C\\. Note that the treatment is identical for the perpendicular component of 


attenuated backscatter: 
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The transmitter-to-receiver overlap function is described by *T(r); however, all regions examined 
in this paper were chosen such that this function is not range dependent and is therefore absorbed 
into the system constants for subsequent equations. 


The accuracies of the CALIOP level 1 products and many of the level 2 products depend 
critically on the accuracy of the calibration of the 532 nm parallel attenuated backscatter profiles. 
The level 1 attenuated backscatter products are used to retrieve the level 2 lidar products, which 
include vertical profiles of aerosol/cloud backscatter and extinction, aerosol/cloud layer base and 
top heights, and integrated aerosol/cloud layer parameters (e.g., aerosol optical depth). Also, 
since the 1064 nm channel is calibrated relative to the calibrated 532 nm signal (Hostetler et al., 
2006, Vaughan et al., 2010), the calibration of all data products are fundamentally dependent 
upon the 532 nm parallel channel calibration. The procedure for calibrating the 532 nm parallel 
attenuated backscatter is described in detail in Powell et al. (2009). Unless otherwise noted, this 
study uses the version 3.01 dataset (V3.01), released in early 2010. 


During nighttime measurements, the CALIOP 532 nm parallel attenuated backscatter is calibrated 
by determining the ratio between the measured signal (i.e., in digitizer counts) at a set altitude to 
the total backscatter estimated for that altitude from an atmospheric model (Powell et al., 2009; 
Hostetler et al., 2006; Russell et al., 1979). The key to accurate calibration is to choose an 
altitude for which the atmospheric backscatter can be accurately estimated and the lidar signal has 



sufficient signal to noise ratio (SNR) and linearity. Solving equation 1 for the parallel calibration 


constant at a given altitude region, r c , yields: 


C„ = 


rlp il( r c) 
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The CALIOP 532 nm parallel channel is calibrated by averaging the signal over the 30-34 km 
range, where aerosol loading is assumed to be low and there is still sufficient molecular 
backscatter to produce a robust signal (Hostetler et al., 2006; Powell et al., 2009). The current 
CALIOP calibration algorithm assumes only molecular backscatter in the calibration region, 
which is computed using molecular number density profiles derived from the NASA Global 
Modeling and Assimilation Office (GMAO) model (Rienecker, 2008) and a known Cabannes 
backscatter cross-section (She, 2001, Hostetler et al., 2006). Estimates of the 532 nm parallel 
channel calibration factor are determined for contiguous 55 km segments of the night side of the 
CALIOP orbit and are smoothed via a 27-segment (1485 km) running average to reduce noise. 
The running average calibration value is used to calibrate the nighttime profiles for the level 1 
attenuated backscatter products. 

For daytime measurements, the high solar background dominates the clear air signal and the 
subsequent low SNR in the 30-34 km region prevents calibration of the parallel channel via the 
molecular normalization technique described above. Instead, the daytime signals are calibrated 
with respect to the calibration established during the previous nighttime orbit. If the magnitudes 
of the nighttime and daytime signals are uniformly proportional with respect to latitude for a 
given target (e.g., cloud-free regions between 8 km and 12 km in altitude), then the nighttime 
calibration coefficient can be readily adapted for daytime conditions via a constant, empirically 
derived scale factor. This approach was used in the initial release of the CALIPSO data, where 



the daytime calibration was estimated by a linear interpolation of the calibration coefficients at 
the endpoints of the two nighttime data segments that bracket each daytime data segment. 
However, it was subsequently shown that this approach led to large errors in the daytime 
calibration due to thermally driven changes in the alignment between the transmitter and the 
receiver as the satellite goes through day-night illumination cycles (Powell et al., 2010). These 
alignment shifts cause the signal levels from identical targets to vary non-linearly over the 
daytime portion of the orbits and thus preclude the use of a constant scale factor to transfer 
calibration from nighttime to daytime measurements. To overcome these instrument anomalies, a 
time-dependent set of scale factors is applied to the nighttime calibration. The scale factors are 
derived using attenuated scattering ratios calculated over cloud-free regions between 8 km and 
12 km in altitude. The goal of this calibration is to produce daytime clean air attenuated 
scattering ratios that are essentially identical to the nighttime clean air attenuated scattering ratios 
at the same altitude and latitude. In the version 2 data releases (V2.0x), the scale factors are 
applied using a five-point piecewise-linear interpolating function (Powell et al., 2008). The 
version 3 data release improves upon this calibration scheme by applying a 34-point latitudinally- 
dependent linear interpolating function (Powell, et al., 2010). In part because of the challenges 
experienced with the daytime calibration and the need for external, independent validation of the 
interpolation function, more HSRL validation flights were conducted during daytime and those 
flights were conducted over as wide a latitudinal and seasonal range as was practical within the 
constraints of the B-200 operational parameters and budget. 

The parallel calibration is transferred to the perpendicular channels using data obtained during a 
polarization calibration operation (Hunt et al., 2009). In this process, a pseudo-polarizer is 
inserted into the optical path, providing equal signal to both channels and thereby allowing the 



electro-optical gain between the parallel and perpendicular channels to be determined, K P gr. The 
532 nm perpendicular calibration constant is then written as: 

C ± = K PGR C U (3) 

This paper focuses on the validation of the 532 nm CALIOP calibration by comparison of the 
total attenuated backscatter ( J3' cauop (r ) ), rewritten from equations 1 and 3: 
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Note that the random uncertainty in Kpgr is small (<1.0%) and systematic uncertainties are 
thought to be much smaller than the random uncertainty (Powell et al., 2009). 

Previous studies have identified systematic biases in the CALIOP attenuated backscatter (Powell 
et al., 2009; Vernier et al., 2009). Sources of systematic uncertainty in the CALIOP calibration 
include uncertainties in the model-derived estimate of backscatter in the calibration region, 
nonlinearity in signal response of the CALIOP detection system, polarization cross talk between 
the parallel and perpendicular 532 nm channels, and errors from radiation-induced noise spikes. 
Based on analysis of pre-launch test data and post launch science data, errors due to detector 
nonlinearity and polarization cross talk are considered to be insignificant (Powell et al., 2009). 
Radiation-induced spikes in the 532 nm signals (e.g., due to high energy protons impinging on 
the detector) can affect calibration by causing errors in background subtraction (i.e., subtraction 
of the digitizer offset and constant background signal from background light) and errors in the 
signal level measured in the calibration range. These spikes are detected and eliminated via an 
adaptive filtering technique (Powell, et al., 2009) prior to the calibration process, and therefore 
should not make a significant contribution to the calibration error. 



Systematic uncertainties in the CALIPSO calibration are discussed by Powell et al. (2009) and 
the CALIOP calibration is expected to have systematic uncertainties of approximately 5%. These 
uncertainties are attributed to uncertainties in the molecular backscatter in the calibration region 
(3%), uncertainties in the aerosol backscatter in the calibration region (4%), or unaccounted 
aerosol, molecular, and ozone transmission between the instrument and the calibration region 
(0.5%). The largest uncertainty is due to assumption of negligible aerosol scattering in the 
calibration region. Indeed, Vernier et al., (2009) estimate that CALIOP is systematically 
calibrated low by ~6% due to unaccounted aerosol scattering up to 35 km in the tropics. Vernier 
et al. (2009) suggests calibration in the relatively aerosol-free region between 36 km and 39 km 
identified in both SAGE and CALIOP data, however this is not implemented in CALIOP data 
processing. This aerosol contribution in the calibration region will be addressed in future 
versions of the CALIOP calibration algorithm. 

3.2 HSRL attenuated backscatter 

The airborne HSRL incorporates 532 nm polarization- sensitive elastic backscatter lidar channels 
similar to those on CALIOP as well as the spectrally-filtered molecular backscatter channel at 
532 nm, from which extinction is derived via the HSRL technique. The relative electrical gains 
and optical attenuations between the channels (the “electro-optic gain ratios”) are determined via 
internal instrument calibrations (Hair et al., 2008). The HSRL 532 nm total attenuated 
backscatter profiles are determined from the internally calibrated profile measurements in a four- 
step process. Lirst, the sum of gain-scaled parallel and perpendicular channel signals is divided 
by the gain-scaled molecular channel signal, yielding the unattenuated scattering ratio, R(r), (the 
attenuation terms in the numerator and denominator cancel in the ratio): 
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Second, this unattenuated scattering ratio profile (R(r)) is multiplied by an estimate of the 
molecular backscatter profile computed from the GMAO-derived molecular density profile and a 
Cabannes-only backscatter cross section (Hair et al., 2008), yielding the total (aerosol plus 
molecular) unattenuated backscatter profile (in km sr' 1 ). 

P( r ) = /?,„ (r) + (3 a (r) = R(r)/3 mModel (r) (7) 

The 532 nm total unattenuated backscatter profile is then used to calibrate the sum of the HSRL 
parallel and perpendicular 532 nm channels to produce the of 532 nm total attenuated backscatter 
profile. 




For the HSRL profiles, ry is chosen to be the highest altitude region in the profile for which the 
transmitter-to-receiver optical overlap factor can be considered to be unity. This is typically 
1.5 km to 2 km below the aircraft (6.5 km to 7 km altitude). The resulting profile is equal to the 
unattenuated backscatter at r/. For altitudes below rj, the P' HSRL (r, r f ) is the product of the total 
backscatter and the two-way transmittance from ry to r. 

The reference altitude range for CALIOP attenuated backscatter is 30-34 km, a calibration 
altitude considerably higher than the altitudes measured by the HSRL. The final step in 
conditioning the HSRL data for assessment of CALIOP calibration is to transfer the reference 
altitude of the attenuated HSRL backscatter to the CALIOP calibration reference altitude of 
30 km. This is done by multiplying the attenuated backscatter profile by an estimate of the 
attenuation between the rf and 30 km. 

30 km 

-2 ja(r')c/r' 

r f 
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The transmittance term includes molecular extinction and ozone absorption losses estimated from 
their number densities in the GMAO model between the HSRL 532 nm total attenuated 
backscatter calibration attitude and that of CALIOP. Background aerosol extinction between the 
two altitudes is unknown and is therefore not included, though this term is estimated to be 
insignificant (see section 4.2). However, clouds or thick aerosols above HSRL can lead to large 
errors in the estimated two-way transmittance. To minimize these errors, the CALIOP level 2 
Vertical Feature Mask (VFM) (Vaughan et al., 2004) is used to detect the presence of clouds and 
aerosols above the aircraft and eliminate high altitude aerosol or cloud-contaminated profiles 
from calibration comparisons. Furthermore, a visual inspection of the level 1 attenuated 
backscatter profiles verified proper VFM identification of high altitude clouds or thick aerosols in 
the region above HSRL for all cases presented in this paper. 

An important aspect of HSRL attenuated backscatter calibration is that it is based on internal 
instrument calibration. The internal calibration enables accurate estimation of total (i.e., aerosol 
plus molecular) backscatter, thereby eliminating calibration errors associated with uncertainties in 
aerosol loading at the calibration altitude. This differs from CALIOP calibration and the 
calibration of backscatter lidars in general which either (a) assume, sometimes incorrectly, that 
the aerosol loading in the calibration region is negligible or (b) requires independent information 
on aerosol loading at the calibration altitude. 

3.3 Sample HSRL underflight from 24 September 2006 

An example HSRL and CALIOP 532 nm total attenuated backscatter comparison from 24 
September 2006 is shown in Figure 2. The time series ‘curtain plots' of the CALIOP (a) and 



HSRL (b) 532 nm attenuated backscatter cover a distance of -270 km along the CALIPSO track 
indicated in the flight track map (c). The enhanced backscattering in the CALIOP measurements 
above 13 km on the southern end of the track (less than 34°N) is due to cirrus clouds. As 
mentioned above, the profiles in which these clouds occur were removed from both the CALIOP 
and HSRL mean 532 nm total attenuated backscatter profiles (d). These mean profiles are the 
basis for the difference calculations: 

difference (r) = 100% • - cal, op O) (1 ^ 

P HSRL (A ^30km ) 

Due to the highly variable nature of clouds and thick aerosols in the planetary boundary layer 
(PBL), the altitude range used in the calculation were restricted to cloud-free regions of the free 
troposphere, from the aircraft to the top of the PBL, where aerosol loading is generally much 
lower and much less variable both spatially and temporally. Altitude regions other than the 
“clean air” free troposphere are addressed on select cases in section 5.4 below. 

In this example, the altitude range selected was in the 4-7 km region and yields a mean difference 
of 2.1% with a standard deviation of 5%. For this case, a less conservative choice on altitude 
range was not found to create significant differences. For instance, using altitude range of 1-7 km 
yields similar results (average difference = 2.7% ± 5%). However, the conservative approach of 
staying well above the PBL is followed for the calibration assessment results presented in this 
paper. 

4. Uncertainty due to the comparison methodology 

Potential systematic biases due to the HSRL validation methodology must be discussed to truly 
assess the CALIOP calibration. This ensures the reported differences of CALIOP 532 nm total 



attenuated backscatter relative to that of HSRL are only due to the CALIOP calibration scheme, 
and not due to the calculation of HSRL attenuated backscatter or the scaling of the HSRL data to 
the CALIOP calibration region discussed in 3.2. Differences due to temporal mismatch of the 
HSRL and CALIOP measurements should average to zero given the large number of cases; 
however this is also verified. 

4.1 Uncertainty due to the HSRL calibration of attenuated backscatter 

Systematic uncertainties between HSRL and CALIOP attenuated backscatter profiles could be 
due to a systematic bias in the HSRL calibration of attenuated backscatter (aerosol or molecular 
backscatter in the 7 km calibration region and the transmission to the calibration region), or the 
scaling of HSRL attenuated backscatter to the CALIOP calibration altitude (30-34 km). 

Similar to CALIOP, HSRL derives an estimate of molecular backscatter using the GMAO 
meteorological model to infer molecular number density and an estimate of the Cabannes cross 
section (Hair et al., 2008). The Cabannes backscatter cross section is a straightforward 
theoretical calculation (She, 2001) and the value used by HSRL is consistent with that of 
CALIOP (within 0.09%). Any uncertainty in the molecular backscatter will be dominated by the 
accuracy of the GMAO model which is used by both CALIOP and HSRL in the data retrievals. 
Hair et al., (2008) estimated the systematic error due to temperature and pressure in the HSRL 
molecular backscatter to be less than ~1%. 

Unlike standard backscatter lidars, the HSRL technique implemented does not rely on a 
calibration from the atmosphere. The internal calibration of the airborne HSRL has been 
thoroughly examined and systematic error in the aerosol backscatter profiles is estimated to be 


less than 2.3% (Hair et al., 2008). 



Atmospheric attenuation between the HSRL (aircraft) altitude and the HSRL attenuated 
backscatter calibration region is due aerosol and molecular scattering and ozone absorption 
(equation 9). The HSRL calculation of attenuated backscatter only accounts for the molecular 
scattering. This is accomplished using the GMAO model and an estimate of the Rayleigh 
scattering coefficient, which amounts to approximately an attenuation of 1.6% for nominal HSRL 
altitudes. Note that the molecular attenuation was not accounted for in Powell et al. (2009), so 
the resulting biases presented there are larger than those reported here. The ozone absorption is 
not accounted for in the HSRL calculation of attenuated backscatter because the ozone absorption 
is estimated to be approximately 0.05% for nominal HSRL altitudes. While HSRL measures the 
aerosol scattering ratio in this region, the aerosol extinction (and hence aerosol transmission) is 
not calculated due to system overlap. The aerosol transmission is therefore neglected in equation 
9; however observing large scattering ratios in this region can be used to quality assure the HSRL 
attenuated backscatter. To estimate the maximum bias due to the neglected aerosol attenuation, 
the aerosol extinction is integrated assuming a scattering ratio of 0.05 (equation 5) and a lidar 
ratio of 50 sr. This amounts to a conservative estimate of less than -0.5% systematic uncertainty 
due to unaccounted aerosol transmission. 

4.2 Uncertainty due to the estimate of the two-way transmittance used to scale airborne HSRL 
attenuated backscatter to the CALIOP reference altitude (30-34 km). 

The two-way molecular and ozone transmission between the HSRL attenuated backscatter 
calibration altitude (-7 km) and the CALIOP calibration altitude (30 km) is typically around 87% 
(Equation 10, T^dec,^ + ozone)- The potential systematic biases which could influence the 
transmission term can be due to ozone, molecules, clouds, and aerosol attenuation. Because the 
aerosol and cloud profiles are unknown, they are not accounted for in the transmission term; only 



the transmission from molecules and ozone is included. Systematic errors in the molecular 
transmission (due to accuracy of the GMAO temperature profile) are thought to be less than 1% 
(Russell et al., 1979). The ozone attenuation term is small (T 2 0z0n e) so any bias induced by this 
correction should be negligible. 

The CALIOP VFM is used to remove any profiles with cloud, aerosol, or stratospheric layers 
above the HSRL measurements; however undetected features will introduce biases into the 
comparisons. To bound the systematic uncertainty introduced by undetected high altitude 
aerosols and clouds, CALIOP level 2 (layer product, version 3.01) cloud optical thicknesses 
(COT) were computed above the HSRL altitude for all 86 flight tracks (Figure 3). The minimum 
CALIOP COT above HSRL provides an estimate of the minimum detectable COT and the 
maximum bias introduced by undetected clouds. Although the HSRL validation flights generally 
targeted cirrus free forecasts, a significant amount of COT over a portion of the HSRL track is 
observed on approximately 25% of the flights. The minimum COT found in this dataset was 
0.0016, but to estimate an upper bound of a COT of 0.0125 is estimated from Figure 3 (slightly 
less than the most frequent COT). Accounting for the multiple scattering enhancement in 
CALIOP with a scaling factor of 0.6 (Winker, 2003), this translates to a COT of 0.075 and 
corresponds to a maximum potential bias due to undetected clouds of 1.5%. This is consistent 
with a value estimated from CPL measurements. McGill et al., (2006) determined that the 
CALIOP minimum detectable backscatter (MDB) (5 km resolution) was 1.7 x 10" km" sr" . 
Assuming the feature is detected at 20 km (multiply MDB by 0.25), a lidar ratio of a mean lidar 
ratio of 25 sr (Winker et al., 2009) (multiply by 25 sr), and layer thickness of 1 km (multiply by 
1 km), the minimum detectable optical thickness is estimated to be 0.01. Now accounting for the 
two way transmittance (multiply by 2) and multiple scattering enhancement (multiply by 0.6), we 



2 

estimate the minimum detectable cloud optical thickness over HSRL to be 1.2% (T c ioud=0.978), 
in agreement with the observation of COT over the HSRL underflights. In summary, thin clouds 
not detected by the CALIOP L2 VFM used for cloud screening could introduce a bias no larger 
than 1.5%. 

Neglecting undetected aerosols may also systematically bias the transmission term. However, the 
optical depth in the upper troposphere and stratosphere is generally low (AOD < 0.005) in the 
absence of any large injections of aerosol into the stratosphere (Jager, 2005). In a relevant 
example, Rogers et al. (2009) found the 532 nm AOD above 6.3 km to be 0.01 during the 
Megacity Initiative: Local and Global Research Observations (MILAGRO) campaign. In order to 
estimate the maximum possible bias introduced by unaccounted aerosol transmission, a constant 
scattering ratio of 1.05 and a lidar ratio of 50 sr is assumed in the 30 km to 8 km region and 
yields a maximum undetected AOD of 0.011 (T aerosoi=0.978). This corresponds to a maximum 
uncertainty of 2.2% due to unaccounted aerosol attenuation. 

4.3 Temporal mismatch between the airborne HSRL and CALIOP observations 
The HSRL and CALIOP mean 532 nm attenuated backscatter profiles in the previous sections 
were averaged in time over the entire (cloud-free) spatially coincident HSRL flight track to 
ensure the best SNR in the CALIOP profile. However, the actual HSRL and CALIOP 
coincidence occurred instantaneously, at a single location and time during each flight. Generally, 
the free troposphere “clean air” is spatially and temporally stable for the time and space scale 
considered in this study (Anderson et al., 2003). The stability is tested by evaluating the clean air 
difference as a function of time from the closest point of approach (CPA) (i.e., where the aircraft 
and satellite are coincident in time). This is accomplished by averaging the HSRL 532 nm total 
attenuated backscatter into 20 minute bins centered about the CPA. The spatially corresponding 



CALIOP data were similarly averaged and the mean differences calculated as functions from the 
CPA. Because each flight has a different mean difference and this hypothesis investigates 
relative changes, the mean difference for the entire flight is subtracted from each 20 minute bin 
for each flight. 

Figure 4 summarizes the relative differences for all flights, with the mean and standard deviation 
of the relative differences plotted for each 20 minute bin as well as the number of points in each 
bin (grey bars). No significant discrepancies (within 1%) are observed as the time from the CPA 
is increased. This indicates that the mean difference for each flight is a good approximation for 
the mean difference at the CPA and the temporal averaging of the entire profiles does not 
influence the analysis in a systematic manner. 

4.4 Summary of systematic uncertainties 

All known potential systematic biases in the validation methodology are presented to ensure any 
reported differences of the CALIOP 532 nm total attenuated backscatter relative to HSRL are 
solely due to the CALIOP calibration scheme. Systematic biases could be introduced by biases in 
the HSRL calibration of attenuated backscatter (aerosol or molecular backscatter in the 7 km 
calibration region and the transmission to the calibration region), or the scaling of HSRL 
attenuated backscatter to the CALIOP calibration altitude (30 km). The temporal mismatch 
between HSRL and CALIOP was not found to introduce a systematic bias. 

All potential systematic uncertainties due to systematic biases are summarized in Table 2. The 
dominate uncertainty source is due to the scaling of the HSRL 532 nm total attenuated 
backscatter from the HSRL calibration altitude (7.5 km) to that of CALIOP (30 km). At most 
these uncertainties can lead to a maximum possible systematic uncertainty of 3.7% (in a root sum 



square (RSS) sense). Because these potential uncertainty sources are not known, they are not 
corrected for in this paper. The intent of Table 2 is to rigorously establish the maximum 
systematic uncertainty from all sources that must be considered for CALIOP validation with an 
airborne lidar. 

5. Results 

5.1 Version 3.01 CALIOP 532 nm total attenuated backscatter comparison with HSRL 
As discussed in section 3.3, relative difference profiles are calculated in the clean air region from 
the mean HSRL and CALIOP 532 nm total attenuated backscatter profiles for each flight. The 
large number of coincident measurements allows a statistical comparison to be performed. The 
resulting differences in the clean air regions identified for all 86 flights are summarized in Figure 
5. Due to differences in the CALIOP calibration approach for the day and night segments of the 
orbit, the calibration differences in Figure 5 are plotted separately based on lighting conditions. 

The nighttime difference distribution is a direct assessment of the CALIOP 532 nm total 
attenuated backscatter calibration procedure and shows a mean difference of +2.7% with a 
standard deviation of 6.3%. The daytime difference is slightly larger, with a mean difference of 
+2.9%, and standard deviation of 20%. The larger spread of the daytime distribution is due to 
noise in the solar background. Indeed, we found that the amount of daytime spread relative to the 
nighttime agrees well with the calculated random error due the background signal of CALIPSO 
for these flights (calculated following Liu et al, 2006). 

In order to assess the latitudinal and seasonal dependence of the version 3.01 dataset, the mean 
difference profiles were averaged vertically and horizontally over cloud-free regions of the free 
troposphere to obtain a single difference point for each HSRL underflight and assigned the mean 



latitude for that underflight. The difference from each flight is shown in Figure 6 as a function of 
latitude and colored by month, similar to results shown in Powell et al., (2009). In this analysis, 
each vertical sample was treated as an independent measurement so the error bars represent the 
standard error of the mean (i.e. the standard deviation divided by the square root of the number of 
points). 

As expected from Figure 5, the mean difference of each flight is generally positive in both day 
and night lighting conditions. The mean of these data points are also positive (2.7% nighttime 
and 2.9% daytime, CALIOP lower), with the standard deviations (2.1% nighttime and 3.9% 
daytime), representing the measured variability of the average calibration difference. 

The nighttime flights show a slight seasonal dependence, with higher differences seen in the 
summertime (June and August) months and lower differences in the spring/winter months. These 
differences are most likely due to stratospheric aerosol in the CALIOP calibration region in the 
summertime months for these latitudes (Hostetler et al., 2006) or aerosol loading in the free 
troposphere (Vernier et al., 2010), and are discussed further in Section 0. The daytime 
differences do not have an obvious latitudinal or seasonal dependence, expected due to 
improvements in the CALIOP version 3.01 daytime calibration via the interpolation function. 
Many of the outliers are from flights where only a short segment of the flight track was 
considered (e.g. screened due to cirrus above the aircraft). The smaller horizontal averaging of 
CALIOP data for these cases results in larger error bars. 

5.2 Version 2.0x CALIOP 532 nm total attenuated backscatter comparison with HSRL 

The version 3.01 CALIOP LI product was released in early 2010, with significant modifications 

made to improve the overall operational code; the modification to the version 3.01 daytime 



calibration algorithm was discussed in section 3.1. An example of the improvement from the 
version 2. Ox to the version 3.01 CALIOP 532 nm total attenuated backscatter daytime 
measurements was observed by HSRL on 16 Oct 2008. The version 2. Ox clean air difference was 
20.2% and the version 3.01 difference was 1.9% for this flight, which is by far the largest 
difference observed in the version 2.0x data product. 

The version 2. Ox clean air differences for all HSRL flights are summarized in Figure 6 (c and d). 
The mean nighttime difference for version 2. Ox is nearly identical (2.2%) to that of version 3.01, 
which is expected because the CALIOP nighttime calibration procedures were essentially 
unchanged in version 3.01. The version 2.0x daytime difference (1.3%) was slightly smaller than 
the difference observed in version 3.01; however this is due to more negative differences in the 
version 2.0x comparison. Indeed, the standard deviation of the version 2. Ox comparison is 5.1% 
in the daytime, which is larger than in corresponding standard deviation in the version 3.01 
comparison (3.9%). Overall, many of the version 2.0x daytime flights with differences larger 
(smaller) than 10% (-10%) in Figure 6 show significant improvement in the CALIOP version 
3.01 calibration, such as the 16 Oct 2008 case highlighted in Figure 7. 

5.3 CALIOP 532 nm attenuated backscatter calibration through the 2009 laser switch 
The CALIOP primary laser was intentionally shut down in Febrary 2009 due to a slow loss of 
pressure in the primary laser canister and CALIOP switched to the backup laser (Hunt et al., 
2009). A dedicated series of nighttime HSRL validation flights was performed to ensure 
calibration consistency across the laser transition, with six flights prior to the laser switch and 
five follow up flights to assess the backup laser. Figure 8 shows the flight tracks and differences 
from these dedicated flights. The differences are nearly identical for the primary laser (0.91%) 



and the backup laser (1.2%) and indicate that the CALIOP 532 nm total attenuated backscatter 
calibration was maintained through the laser switch. 

5.4 Altitude dependence of the CALIOP 532 nm total attenuated backscatter? 

To assess whether the CALIOP 532 nm total attenuated backscatter has any altitude dependence 
with respect to that of HSRL HSRL, difference profiles were computed as a function of altitude 
over the full vertical range of the HSRL profiles. While Anderson et al. (2003) indicate that there 
should be little expected variation in aerosol concentration within a few hours on average, the 
largest variations are expected to occur inside the PBL (rather than in the free troposphere) due to 
local sources of aerosol, higher relative humidity, changing meteorological conditions, and other 
similar factors. In order to examine the temporal variability further, the HSRL attenuated 
backscatter profiles were examined on flights that flew a CALIPSO track and then back tracked 
the same path on the return to base. The overlapping tracks were matched up in a gridded 
latitude/longitude space, allowing a direct measurement of how much the attenuated backscatter 
profile changed as a function of time (at multiple matching locations along the track). This 
analysis suggested that the HSRL attenuated backscatters were well correlated (r > 0.9) in the 
PBL with up to at least 45 minutes separation. 

Therefore, a subset of the 86 flight comparison dataset was considered with tighter constraints on 
temporal separation of HSRL and CALIOP. In this subset, only HSRL data within 45 minutes of 
the CALIOP closest point of approach were considered, which is slightly larger than the 30 
minute criteria used by Mona et al., (2009). Furthermore, because clouds are highly variable in 
the atmosphere both HSRL and CALIOP data were screened for clouds at all altitudes using both 
the CALIOP VFM and HSRL cloud detection routine (detected as sharp gradients in the raw 
signals using a three point Haar wavelet covariance transform (Gamage and Hagelberg, 1993)), as 



well as a manual inspection to ensure that only cloud-free profiles were used in the comparison. 

A total of twenty two cases (eight nighttime and fourteen daytime) were found to meet this 
selection criteria. The HSRL and CALIOP mean 532 nm total attenuated backscatter profiles for 
these cases are presented in Figure 9. 

The mean difference profiles from all 22 flights is shown in Figure 10, separated by day and 
night. For plotting purposes, the differences were binned into 500 m bins, with the error bars 
representing the standard error of the mean for each bin. Further averaging the day and night 
profiles, a near-constant difference of 3.1% is observed throughout the profile, which is 
consistent with the “clear air” difference of these 22 flights (2.4%). To estimate the change over 
altitude, a linear regression was performed treating the altitude as the independent variable. The 
slope of 0.011 % per kilometer suggests that over a 10 km range the difference will vary by less 
than 0.1% due to vertical variability and the vertical range selected for each flight does not 
influence the results in cloud-free regions. 

Figure 11 shows the mean CALIOP 532 nm total attenuated backscatter plotted against that of 
HSRL for all altitudes of all 22 cases. Clearly, they are well correlated (r = 0.88) and have a 
slope near unity (slope = 1.03) across a large range of attenuated backscatter values. This 
demonstrates the consistency accuracy of the CALIOP attenuated backscatter over the large 
dynamic range of aerosol loading present in these 22 cases. Note that the lowest attenuated 
backscatter reported by HSRL (-0.6x10" km" sr" ) is limited by the altitude of the aircraft. . 

6. Discussion and Conclusions 

The results and methodology of the most extensive systematic study of the CALIPSO 532 nm 


total attenuated backscatter product to date are presented using an internally calibrated airborne 



HSRL, providing an independent verification of the CALIOP calibration. The version 3.01 
CALIOP 532 nm total attenuated backscatter was found on average to be only slightly (2.7% ± 
2.1%) lower that the HSRL 532 nm attenuated backscatter in clean air during nighttime 
measurements, directly validating the CALIOP calibration algorithm. A slight seasonal 
dependence is observed in the nighttime differences, which is attributed to the stratospheric 
aerosol influence on the CALIOP calibration in the summertime months for the latitudes of 
nighttime HSRL operation. Additionally, this study also found no change in the nighttime 
calibration of 532 nm total attenuated backscatter through the laser change out in early 2009. 
This study also quantitatively assesses the daytime version 3.01 CALIOP 532 nm total attenuated 
backscatter coefficients, which are found on average to be only slightly (2.9% ± 3.9%) lower 
than the HSRL 532 nm attenuated backscatter measurements. Furthermore, the daytime 
calibration values exhibited no obvious latitudinal or seasonal dependence (within -10%), 
indicating accurate performance of the version 3.01 daytime interpolation function. 

The comparison methodology presented here is thoroughly examined and intended to establish a 
model for satellite lidar validation with an airborne lidar. Systematic biases considered in this 
methodology include the following: HSRL calibration of attenuated backscatter due to errors in 
the aerosol and molecular backscatter in the calibration region and transmission to the HSRL 
calibration region; biases due to scaling the HSRL attenuated backscatter to the CALIOP 
calibration region; biases due temporal and spatial offsets of the HSRL and CALIOP 
measurements. Overall, these errors amount to a systematic uncertainty no larger than 3.7%. 
This implies that the CALIOP 532 nm total attenuated backscatter profiles are well calibrated 
within the ability of HSRL to measure. 



In addition to aircraft studies previously discussed (McGill et al., 2006), ground based lidars have 
also offered an assessment of CALIOP attenuated backscatter profiles. Kim et al. (2009) 
qualitatively shows good agreement between an elastic backscatter lidar and CALIOP attenuated 
backscatter. In another study, Mamouri et al., 2009 used a ground based Raman lidar for 
nighttime validation and an elastic backscatter lidar for daytime validation of the CALIOP 532 
nm total attenuated backscatter. Note that a Raman lidar has an advantage for validation over an 
elastic backscatter lidar because the Raman lidar technique, like the HSRL technique, provides a 
more direct measurement of the extinction profile. However, comparisons of nadir (CALIOP) 
and zenith (ground-based) lidar profiles can have higher uncertainties due to differences in 
viewing geometry. For example, Ansmann (2006) found that Klett retrievals from zenith- 
viewing lidar could differ in extinction and backscatter profiles by as much as 20% from a nadir- 
viewing lidar observing the same scene due to the viewing geometry. 

Two ground based Raman lidar studies have provided quantitative assessment of the CALIOP 
532 nm attenuated backscatter calibration and show promising results. In clear sky conditions, 
Mona et al. (2009), found CALIOP to bias slightly low in the free troposphere and very low in 
the PBL (below 2.5 km). Pappalardo et al., (2010) show similar results from comparisons 
between several EARLINET lidars and CALIOP with similar results. Additionally, Powell et al. 
(2009) present the first nighttime HSRL validation flights through 2008, which are a subset of 
data presented here. However, the results from Powell at al. (2009) did not account for molecular 
transmission in the HSRL calibration (see section 4.1) which resulted in a systematic high bias in 
HSRL by approximately 1.6%. Table 3 summarizes the published lidar studies which 
quantitatively assess CALIOP level 1 attenuated backscatter to date. 



The differences found in this study agree in the free troposphere with the results from previous 
studies. However, this study found good agreement between HSRL and CALIOP 532 nm total 
attenuated backscatter inside the PBL (below ~3km) while previous studies suggest that CALIOP 
may be biased low in this region. In section 5.4 the vertical profile analysis clearly demonstrates 
that the CALIOP attenuated backscatter profiles are linear through the vertical range investigated 
and the assessment in the free troposphere applies to the entire vertical profile. Note that these 
previous quantitative studies were conducted using ground based lidars. Due to differences in the 
nadir (CALIOP) versus zenith (ground-based) viewing geometries and spatial mismatch, such 
comparisons may lead to ambiguous results. Discrepancies between nadir- and zenith- viewing 
approachs may be due to spatial offset of the ground station and the CALIPSO track or temporal 
averaging of the CALIOP data around the overpass, coupled with influences of local sources and 
complex terrain between the satellite track and the ground station (Mona et al., 2009). The 
largest discrepancies due to these considerations are expected appear in the PBL. Indeed, 
consider the average minimum distance presented by Mona et al., (2009) was 66.5 km. 
Assuming a uniform 20 km/hr wind (also assuming directly between the CALIPSO track and the 
lidar), the minimum equivalent temporal separation between the measurements is over 3 hrs, 
when the aerosol correlation starts to decrease (Anderson et al., 2003). Another potential error 
source is the estimation of a nadir viewing attenuated backscatter profile from a zenith viewing 
lidar (Pappalardo et al., 2010), with potential errors at the lower altitudes due to inaccuracies in 
the lidar overlap function (Wandinger and Ansmann, 2002). We stress that the airborne HSRL 
measurements presented here are in agreement with previous (ground based) studies in the free 
troposphere where ground based overlap and local aerosol influences are less of an issue than in 
the PBL. We therefore conclude that the there is no altitude dependence of the calibration. 



This study further expands the scope of previous validation efforts by including both night and 


day comparisons that cover a wide range of latitudes, thus greatly increasing the number of 
validation cases of CALIOP 532 nm attenuated backscatter with considerably smaller systematic 
uncertainties. Finally, this study takes advantage of the recent release of the CALIOP version 
3.01 data, and provides validation of both version 2. Ox and version 3.01 products. Looking at the 
differences in the CALIOP version 2. Ox and version 3.01 532 nm total attenuated backscatter, the 
mean version 3.01 and version 2. Ox nighttime differences were found to have similar values and 
standard deviation. The daytime version 3.01 differences shows smaller discrepancies compared 
to HSRL than that of version 2. Ox. 

The HSRL has acquired a substantial CALIOP validation dataset over a large temporal and 
spatial range. Nadir-viewing observations made along the CALIPSO ground track via the self- 
calibrated HSRL technique provide, to our knowledge, the best means by which to assess 
CALIOP calibration accuracy and the level products. This rich dataset will be used in the future 
to continue to assess the CALIOP level 1 532 nm and 1064 nm products as well as the level 2 
products, which will be the subject of a future study. 
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Figure 1: HSRL flight tracks of CALIPSO underpasses flown from June 2006 through 2010 

with approximate mission location noted. Blue flight tracks indicates night time flights. 

Figure 2: The 532 nm attenuated backscatter curtains referenced to 30 km from (a) CALIPSO (5 
pt temporal running average) and (b) HSRL on 24 September 2006 from 08:00 to 09:15 UTC. 
The white vertical line indicates the closest point of approach (CPA), 08:20 UTC, and the 
magenta horizontal line indicates the B-200 flight altitude. The flight track map (c) shows the 
HSRL flight path (blue) and CALIOP (red). Mean HSRL (blue) and CALIOP (red) 532 nm total 
attenuated backscatter are shown in (d). 

Figure 3: Histogram of cloud optical thicknesses for clouds above 7.5 km and consequently 
above the maximum height observed by HSRL from all flights. The vertical black line 
corresponds to a COT of 0.0125. 

Figure 4: The relative difference as a function of time from the CPA. Relative difference were 
binned into 20 minute bins with the mean and standard deviation of each bin reported (black) and 
the number of points in each bin (grey bars). 

Figure 5: Distribution of the average 532 nm attenuated backscatter difference in the clean air 
regions identified for all night (blue) and day (red) flight. 

Figure 6: Version 3.01 532 nm total attenuated backscatter differences for (a) night and (b) 
daytime lighting conditions as functions of latitude and month (color). Each data point represents 
the average clean air difference over an entire HSRL flight. Version 2.0x differences (discussed 
in 5.2 below) are also shown for night (c) and day (d) 

Figure 7: Time series (a and b) same as Figure 1 except for daytime flight on 16 Oct 2008 over 
North Carolina and Virginia. Average 532 nm total attenuated backscatter profiles (cloud 
screened) for CALIOP version 2. Ox and version 3.01 are also shown compared to HSRL (c and 
d). 

Figure 8: Flight track ap showing locations of HSRL underflights of CALIPSO between January 
and April 2009 (a) and the corresponding 532 nm attenuated backscatter differences. Night flights 
show the consistent calibration of CALIOP primary (black) and backup (blue) laser. 

Figure 9: Mean CALIOP (red) and HSRL (blue) 532 nm attenuated backscatter for all 22 cases 
meeting the criteria in the text for an assessment of the complete CALIOP vertical profile. 

Figure 10: The mean difference between HSRL and CALIOP for all 22 cases presented in Figure 
9, separated by day (red) and night (blue). The black line represents the linear regression of all 
22 cases (treating altitude as the dependent variable). 

Figure 11: Scatterplot of HSRL and CALIOP 532 nm total attenuated backscatter values from 
the 22 cases presented in Figure 9. The white dashed line indicates the one to one line and the 
solid white line is a bilinear fit. The colorbar indicates the number of points in each bin. 



Mission 

Date Range 

Number of 
CALIOP flights 

Number of hours 
on CALIOP track 

CC-VEX 

14 June 2006 - 17 Aug 2006 

11 

16.2 

TexAQS -GOMACCS 

28 Aug 2006 - 28 Sept 2006 

10 

13.8 

CHAPS 

03 June 2007 - 26 June 2007 

8 

10.9 

CATZ 

19 July 2007- 11 Aug 2007 

4 

7.6 

Caribbean 

24 Jan 2008 - 03 Feb 2008 

7 

13.2 

ARCTAS (spring) 

01 April 2008 - 19 April 2008 

12 

17.5 

ARCTAS (summer) 

14 June 2008 - 10 July 2008 

11 

10.3 

Nighttime Calibration 

22 Jan 2009 - 17 April 2009 

11 

15.9 

RACORO 

17 June 2009 - 26 June 2009 

3 

4.0 

Other 

2007-2009 

9 

6.3 

Total 


86 flights 

115.7 hrs 


Uncertainty in the HSRL calculation of attenuated backscatter 


HSRL Molecular Backscatter 

±1% 

HSRL Aerosol Backscatter 

±2.3% 

Transmission to calibration region 

±0.5% 

Uncertainty in scaling the HSRL attenuated backscatter to 30 km 


Molecular Scattering 

~0% 

Ozone Absorption 

~0% 

Aerosol Scattering 

0-2.2% 

Cloud Extinction 

0-1.5% 

Maximum RSS Uncertainty 

3.7% 







Study 

Number of 
cirrus free cases 

532 nm attenuated 

backscatter difference 
[%] 

Altitude range 

Pappalardo et al. (2010) 

46, night 

4.6% ± 50% 

1-10 km 

Mona et al. (2009) 

11, night 

24+ 20% 

<2.5 km 



2 ± 12% 

3-8 km 

Mamouri et al. (2009) 

12, day 

10+12% 

3- 10km 



34+34% 

l-3km 


15, night 

4+6% 

3- 10km 



15+16% 

l-3km 

Powell et al. (2009) 

9, night 

5% 

~3 km - 7 km* 

This study (version 2) 

56, day 

1.3 + 5.1% 

~3 km - 7 km* 


20, night 

2.2 + 1.7% 


This study (version 3) 

66, day 

2.9 + 3.9% 

~3 km - 7 km* 


20, night 

2.7 + 2.1% 



Table 1: Summary of flights and hours along the CALIPSO track for the field missions 

containing CALIOP validation components 

Table 2: Summary of systematic uncertainty sources that may influence the validation of the 
CALIOP 532 nm total attenuated backscatter with the airborne HSRL measurements 


Table 3: Summary of lidar studies to date quantitatively validating the CALIOP level 1 532 total 
attenuated backscatter product. Note that all differences are reported in terms of ‘Lidar - 
CALIOP’ and that all previous studies used the version 2 data release. 

*altitude ranges varied, but no dependence noted. 
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